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Abstract 



We have performed shell-model calculations for the N = 126 isotones 210 Po, 
211 At, and 212 Rn using a realistic effective interaction derived from the Bonn 
A nucleon-nucleon potential by means of a G-matrix folded-diagram method. 
The calculated binding energies, energy spectra and electromagnetic prop- 
erties show remarkably good agreement with the experimental data. The 
results of this paper complement those of our previous study on neutron hole 
Pb isotopes, confirming that realistic effective interactions are now able to 
reproduce with quantitative accuracy the spectroscopic properties of complex 
nuclei. 

PACS number(s): 21.60.Cs; 21.30.Fe; 27.80. +w 



1 



Typeset using REVTgX 



I. INTRODUCTION 



During the past few years, we have studied a number of nuclei around doubly magic 100 Sn 
and 132 Sn |l]-[5] in terms of the shell model employing realistic effective interactions derived 
from the meson-theoretic Bonn A nucleon-nucleon (NN) potential ||. In these studies we 
have considered nuclei with few valence particles or holes, their properties being of special 
interest for a stringent test of the basic ingredients of shell-model calculations. The aim of 
our work is to assess the ability of realistic effective interactions to provide a quantitative 
description of nuclear structure properties. This is in fact a key point to understand if the 
time has come to make the shell model a truly microscopic theory of nuclear structure. 

As is well known, the first step in this direction was taken more than thirty years ago by 
Kuo and Brown |7|] who derived an sd-shell effective interaction from the Hamada- Johnston 
potential || . Later on, an effective interaction for the lead region was derived |J by Kuo and 
Herling (KH) from the same potential. Since then, however, substantial progress has been 
made in both the development of high-quality NN potentials and the many-body methods 
for calculating the matrix elements of the effective interaction. 

As regards the first point, modern potentials reproduce quite accurately all the known 
NN scattering data. A review of recent developments in the field of NN potentials is given 
in Refs. [|TU,rT|. We only recall here that two potentials which fit equally well the NN 



data up to the inelastic threshold may differ substantially in their off-shell behavior. Thus, 
different NN potentials may produce somewhat different nuclear structure results. 

As for the second point, an accurate calculation of the Brueckner G matrix is now feasible 
while the so-called folded-diagram expansion yields a rigorous expression for the model-space 
effective interaction V e R. The main aspects of the above derivation of V e s are described in 



Refs. plpl 



Based on these improvements, a new generation of realistic effective interactions has 
become available, fostering renewed interest in realistic shell-model calculations. It is in 
this context that our recent studies of medium-mass nuclei are framed. The remarkably 
good agreement between theory and experiment obtained for these nuclei has challenged us 
to perform the same kind of realistic shell-model calculations for heavy-mass nuclei. In a 
previous work [|14| we considered the neutron hole isotopes 206 > 205 > 204 pb. Here, we present 
the results of a companion study of the iV = 126 isotones, focusing attention on 210 Po, 
211 At, and 212 Rn. These nuclei, with two to four protons in the Z = 82 — 126 shell, offer the 
opportunity to further test our realistic effective interactions in the lead region. 

The N = 126 nuclei, as well as the lead isotopes, have long been the subject of both 
experimental and theoretical studies. From the experimental point of view, these stable or 
near-stable nuclei have been extensively investigated and a rather large amount of exper- 
imental data is available for them. On the other hand, the good doubly magic character 
of 208 Pb has motivated many shell-model calculations in this region. In all the calcula- 



tions performed so far, however, phenomenological interactions have been used [15,16], the 



only notable exception being the pioneering work of McGrory and Kuo |T7| , where the KH 
interaction was employed. 

The outline of the paper is as follows. In Sec. II we give an outline of our calculations, 
including a brief review of the derivation of the effective interaction. In Sec. Ill we present 
the results obtained for binding energies, energy spectra and electromagnetic properties, 
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comparing them with the experimental data. Sec. IV contains a discussion and a summary 
of our conclusions. 



II. OUTLINE OF CALCULATIONS 

We assume that 208 Pb is a closed core and let the valence protons occupy the six single- 
particle (sp) orbits O/tg/2, I/7/2, I/5/2, 2p 3 / 2 , and 2pi/ 2 . We take the sp energies from 
the experimental spectrum of 209 Bi Hl8| . They are (in MeV): eh g/2 = 0.0, e/ 7/2 = 0.896, 
e il3/2 = 1.609, e h/2 = 2.826, e P3/2 = 3.119, e Pl/2 = 3.633. 

As already mentioned in the Introduction, we use in the present calculation a realistic 
effective interaction derived from the Bonn A free NN potential. Let us now outline our 
derivation of V e g. Because of the strong repulsive core contained in the Bonn A potential, 
which is a feature common to all modern NN potentials, the model-space G matrix corre- 
sponding to the chosen Vnn must be calculated first. The G matrix is defined |L9| by the 
integral equation: 

G{u>) = V + VQ 2 ±-—Q 2 G(u), (1) 

LO — C^ 2 i V2 

where V represents the NN potential, T denotes the two-nucleon kinetic energy, and u is an 
energy variable (the so-called starting energy). The two-body Pauli exclusion operator Q 2 
prevents double counting, namely the intermediate states allowed for G must be outside of 
the chosen model space. Thus the Pauli operator Q 2 is dependent on the model space, and 
so is the G matrix. The operator Q 2 is specified, as discussed in Ref. [ID], by a set of three 



numbers (77.1,712,713) each representing a shell-model orbital (we number the orbits starting 
from the bottom of the oscillator well; for instance, the orbit 0^5/2 is denoted as orbit 4 
and O/ixi/2 as orbit 16). Note that in Eq.(l) the Pauli exclusion operator Q 2 is defined in 
terms of harmonic oscillator wave functions while plane-wave functions are employed for the 
intermediate states of the G matrix. 

Since the valence-proton and -neutron orbits outside 208 Pb are different, our Q 2 operators 
for protons and for neutrons are different and consequently our G-matrix calculation is 
considerably more complicated than in the case when the two operators are the same. In 
the present calculation we have fixed (m, n 2 , n 3 ) = (22,45,78) for the neutron orbits, and 
(ni, n 2 , ris) = (16, 36, 78) for the proton orbits. Our procedure for calculating the G matrix is 
outlined below. We first calculate the free G matrix Gf in a proton- neutron representation, 
Gp being defined by 

G F = V + V—Gp, (2) 

e 

with e = (u — T). Note that Gf does not contain the Pauli exclusion operator and hence 



its calculation is relatively convenient. Then we calculate the Pauli correction term |19| , |20 



where P 2 = 1 — Q 2 , separately for protons and for neutrons. Finally the full G matrix as 
defined by Eq.(l) is obtained as 



3 



G = G F + AG. 



(4) 



For the harmonic oscillator parameter huo we use the value 6.88 MeV, as obtained from 
the expression htu = 45A~ 1//3 — 25A~ 2 ^ 3 for A = 208. Note that in the present work we 
have chosen the value of n 2 so as to include three harmonic oscillator shells above the Fermi 
level denoted by n\. In earlier works on light- and medium-mass nuclei, n 2 was fixed by 
taking into account only two major shells above the ni-th orbit. For instance, a common 
choice for sci-shell calculations is n 2 =10 with rii=3 For the N = 126 isotones, however, 
the present choice is more appropriate. In fact, these calculations as well as those for lead 
isotopes jl4]2~l]l have led to a substantially better agreement with experiment when in the 
derivation of the effective interaction n 2 has been increased from two to three shells above 
the ni-th orbit. 

Using the above G matrix we then calculate the irreducible vertex function Q-box, which 
is composed of irreducible valence-linked G-matrix diagrams through second order in G. 
These are precisely the seven first- and second-order diagrams considered by Shurpin et al. 



22]. The effective interaction can be written in operator form as 



V. 



eff 



Q-Q' Q + Q' Q Q-Q' Q Q Q + 



(5) 



where Q is the Q-box, and the integral sign represents a generalized folding operation |23 
Q' is obtained from Q by removing terms of first order in the reaction matrix G. After 
the Q-box is calculated, the energy-independent V e e is then obtained by summing up the 
folded-diagram series of Eq.(5) to all orders using the Lee-Suzuki iteration method ||24j| . This 
last step can be performed in an essentially exact way for a given Q-box. 

Once the effective interaction has been derived, the shell-model calculations are carried 



out employing the OXBASH code [25 



As regards the electromagnetic observables, we have calculated them by making use of ef- 
fective operators 
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which take into account core-polarization effects. More precisely, by 



using a diagrammatic description as in Ref. pq| , we have only included first-order diagrams 
in G. This implies that folded-diagram renormalizations are not necessary ||271| . 

As is well known, the nuclear magnetic properties may be significantly affected by mesonic 
exchange currents. An estimate of their contribution for nuclei in the vicinity of 208 Pb has 



been given in Refs. [28,29]- This amounts to renormalizing the gyromagnetic factor gi from 
the bare value of gi = 1 to 1.155 and g s from 5.586 to 5.699. We have made use of these 
values in our calculation of the effective Ml operator. 



III. RESULTS AND COMPARISON WITH EXPERIMENT 

In this section the results of our calculations for the three nuclei 210 Po, 211 At, and 212 Rn 
are presented and compared with experiment. The energy spectra are presented separately 
for each isotone in the three following subsections. For 210 Po a detailed comparison between 
calculated and observed spectroscopic factors is also reported. The last subsection is devoted 
to the discussion of the electromagnetic properties of all three nuclei. 

The calculated ground-state binding energies relative to 208 Pb are compared with the 
observed values |Kj in Table I. The mass excess value for 209 Bi needed for absolute scaling 
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of the sp levels was taken from [ 30 1 . As regards the Coulomb interaction between the 



valence protons, we have assumed that it gives a contribution proportional to the number of 
interacting proton pairs, namely Eq = Vp (n is the number of valence protons). The 

strength Vq has been taken to be 270 keV, which is the value of the matrix element of the 
Coulomb force between two hg/2 protons with J = 0. 

From Table I we see that a very good agreement with experiment is obtained for all three 
nuclei. In fact, the calculated binding energy for 210 Po falls practically within the error bar, 
while the other two calculated values differ by less than 100 keV from the experimental ones. 



A. Spectrum of 210 Po 



The experimental [jTS] and theoretical spectra of 210 Po are compared in Fig. 1. Here 
all the calculated and experimental levels up to 3.3 MeV are reported, while in the higher- 
energy region the negative-parity states have been excluded (in the experimental spectrum 
we have also omitted three states without spin and parity assignment). The calculated 
states which are shown in Fig. 1 are all the 44 states arising from the configurations hg, 2 , 
hg/2 f 7/2, h 9/2 ii 3 /2, fy 2 , «i3/ 2 > ^9/2/5/2, and h 9/2 p 3 /2- In the energy region 3.5-4.5 MeV we find 
the eight negative-parity states of the 77/2*13/2 configuration, but for them, as mentioned 
above, we have not tried to establish any correspondence with observed levels. In fact, in 
this energy interval negative-parity states have been observed which cannot be described 
within our model space. As an example, we mention the 12~ and 13~ states, which cannot 
be constructed in our model space, and the three 5~ levels observed at 3.43, 3.70, and 3.71 
MeV to which corresponds only one calculated 5~ state at 3.85 MeV. These experimental 
"extra" levels arise from core excitations, and in some cases significant admixtures of these 
excitations and two-particle model-space states are likely to occur. 

In Fig. 1 we see that the calculated spectrum is characterized by four groups of levels: 
the first one up to 1.6 MeV, the second between 2.2 and 2.5 MeV, the third between 2.7 and 
3.3 MeV, and the fourth above 3.8 MeV. The five levels of the first group are dominated 
by the h\i 2 configuration, while the second group contains all the members of the /19/2/7/2 
multiplet, besides the + state arising from the fy 2 configuration. The other three states of 
this latter configuration together with all the states arising from the ^.9/2213/2 configuration 
are in the third group. All the states in these three groups, with few exceptions, are almost 
pure, the percentage of the dominant configuration ranging from 95 to 100%. Only for the 
ground state and the 0^ and 2% states the contribution coming from configurations other 
than the dominant one is particularly significant, the percentage of such configurations being 
21, 28, and 16%, respectively. The 17 levels arising from the configurations if 3 / 2 , ^9/2/5/2, 
and I19/2P3/2 all he in the energy interval 3.8-4.8 MeV. It should be mentioned, however, that 
for most of these states the wave functions are not quite pure. In particular, we find that a 
significant admixture of the three above configurations is present in the even J states. 

Up to 3.3 MeV excitation energy each state of a given J n in the calculated spectrum 
can be unambiguously associated with an observed level, the only exception being the 2g~ 
state at 2.95 MeV excitation energy. However, two levels with no angular momentum and 



parity assignment have been observed at 2.66 and 2.87 MeV, and in Ref. |3T| it is suggested 
that the 2.87-keV 7 ray measured in the 209 Bi(t, 2n) 210 Po reaction is a good candidate for 
the 2% — > 0+ transition. The experimental 37 state at 2.39 MeV, as well as the the 57 
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and A 2 states at 2.91 and 3.11 MeV, respectively, have no theoretical counterpart. In fact, 
the first one reflects the collective nature of the octupole 3~ state at 2.61 MeV in 208 Pb, 
while the other two levels arise from the neutron particle-hole configuration v(g 9 / 2 Pi/ 2 ) P^l > 
and therefore cannot be described within our model space. It should be noted that each 
of the two above and 4 2 levels lies very close in energy to the state with the same J n 
originating from the ^9/2*13/2 configuration. Therefore it cannot be excluded, as we shall see 
when discussing the spectroscopic factors, that some mixing occurs between single-particle 
and core-excited states. 

Above 3.8 MeV only 10 out of the 17 levels arising from the configurations i1 3 / 2 , ^9/2/5/2, 
and h 9 / 2 p 3 / 2 have been experimentally identified. For all of them, except the (7 + ,4 + ) state 
at 4.55 MeV, a correspondence with states predicted by the theory can be safely established. 
As for the (7 + , 4 + ) state, it may be associated with either the 7 2 or the 4f calculated states, 



which lie at 4.48 and 4.52 MeV excitation energy, respectively. In a very recent work |32 
the assignment 7 + has been proposed for the experimental level at 4.55 MeV and a new level 
with J 77 = (4+) has been identified at 4.54 MeV. 

As regards the quantitative agreement, we see from Fig. 1 that it is very satisfactory, 
the discrepancies between calculated and experimental excitation energies being less than 
100 keV for most of the states. More precisely, including the level at 2.87 MeV (identified 



as a J 71 " = 2 + state) as well as the J n = 4 + and 7 + of Ref. [32], 37 observed levels have been 
associated with states predicted by the theory, and only for seven of them the experimental 
and calculated excitation energies differ by more than 100 keV. The rms deviation o |33| 
relative to these 37 levels is 87 keV. 



In Ref. the 209 Bi( 3 He,rf) 210 Po and 209 Bi( 4 He,t) 210 Po reactions have been studied and 
the single-proton strengths of transitions to various excited levels in 210 Po have been ex- 
tracted from the measured cross sections. These observed strengths are compared with the 
calculated values in Table II, where we also list the experimental and theoretical excitation 
energies. The experimental uncertainties appearing in Table II are only statistical and the 
numbers in parentheses, which correspond to levels not fully resolved, were extracted by a 
peak fitting procedure (see Ref. The theoretical spectroscopic factor S is defined as 

s£ p = 2jVtK 310r >> r ' P 11 ^ irBi ' J ? = 9/2 "' 9s) ^ 

where we assume that the ground state of 209 Bi is a single hg/ 2 proton outside the doubly 
magic 208 Pb. The label (3 specifies states of 210 Po with the same J 71 . 

From Table II we see that for almost all the states of the three low-lying multiplets the 
agreement between theory and experiment is very good. Actually, a significant discrepancy 
exists only for the 11~ and the second 5~ state. However, it was suggested in Ref. |34J that 
the level at 2.85 MeV was an unresolved doublet with J n = 11~ and 3~, as it has been found 
later to be the case [|L8|]. Thus the observed strength, 3.25, attributed to the 11~ state has 
to be compared to the sum of the two calculated strengths relative to the 11~ and 3~ states, 
which are 2.30 and 0.69, respectively. As for the 5~ state, our calculation overestimates the 
experimental value. Part of the single particle strength is contained in the first 5~ state, 
which, as mentioned above, is not predicted by the theory, being essentially a core-excited 
state. 

In the region above 3.8 MeV the values of the measured strengths (1 = 1 and 3) are 
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generally smaller than those relative to the states of the three low-lying multiplets. On the 
other hand, as it was already pointed out at the beginning of this section, the calculated 
wave functions of several states in this region show a strong admixture of the configurations 
^13/2 > ^9/2/5/2, and h 9 / 2 p 3 / 2 . Thus, a comparison between theory and experiment may provide 
a test of the calculated percentages of the ^9/2/5/2 and /19/2P3/2 configurations (obviously, 
the contribution of the i\ 3 / 2 configuration is not determined directly from the measured 
strengths). It should also be noted that the experimental data do not allow to distinguish 
between p 3 / 2 and py 2 transfers. We have found, however, that a small component of the 
/19/2P1/2 configuration is present only in the 4^ and 4$ states. From Table II we see to the 
observed strengths of the states above 3.8 MeV are quite well reproduced by the theory. Note 
that the level at 4.55 MeV excited via f'5/2 transfer is likely to correspond to an unresolved 
doublet with J n = 4 + and 7 + (see discussion above). In this case the measured strength, 
1.83, should be interpreted as the sum of the calculated strengths 1.50 and 0.20 relative to 
the 7 2 and 4g~ states, respectively. In this connection, it should be pointed out that the 
observed strength of the level at 4.55 MeV excited via I = 1 transfer and assigned J 71 " = 4 + 



34| is also well reproduced by our calculation. 



B. Spectrum of 211 At 

The experimental flTj| and theoretical spectra of 211 At are compared in Fig. 2, where 



all the observed levels up to 3.3 MeV excitation energy are reported. In the calculated 
spectrum all the levels up to about 2.0 MeV are included while in the higher-energy region 
only the states which can be associated to the observed ones are reported. For the sake of 
completeness all the calculated excitation energies up to 2.7 MeV are listed in Table III. 

From Fig. 2 we see that a one-to-one correspondence can be established between the 
experimental and calculated levels up to 1.5 MeV, the only exception being the experimental 
(|, y , y) level at 1.35 MeV which can be associated to either the (| ) 2 or (y + )i calculated 
state. As regards the two observed levels with no firm spin assignment at 1.12 and 1.23 
MeV, we propose the assignment J 71 = y and y , respectively. 

Above 1.5 MeV many more levels than the observed ones are predicted by our calcu- 
lations. In particular, in the energy interval 1.5-2.0 MeV we find 12 levels, only three of 
them having the experimental counterpart. These three states have J n = (|) _ , (y ), and 
(|)~ and can be identified with the calculated states with the same angular momentum and 
parity at 1.82, 1.97, and 2.04 MeV, respectively. It should be noted that we predict the 
existence of a lower- lying | state at 1.85 MeV. 

As regards the states above 2.0 MeV, we identify the experimental J n = | , § level at 
2.06 MeV with the calculated one with J n = | at 2.11 MeV, thus confirming the tentative 
assignment of Ref. |35|]. As far as the lowest | + state is concerned, the calculated energy 
is 2.78 MeV, namely 300 keV higher than that of the observed one. However, this state, 
which was populated in a first-forbidden [3 + decay of the ground state of 211 Rn P5| , has 



been interpreted as a core-excited 212 Rn (^(vrs!^) -1 state and, as such, is not expected to 
be adequately reproduced within our model space. In Ref. the nature and assignment 
of the J = |, § level at 2.65 MeV was also discussed and it was suggested that it originates 
from a mixing of the two 212 Rn ®{iid 3 / 2 )~ 1 and 211 At(gs) ® 210 Pb(3~) core-excited states. 
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We may only mention here that up to 3.5 MeV our calculations predict the existence of the 
three ~ states reported in Fig. 2 and of no other ^ aside the above mentioned one. 

The quantitative agreement between our results and experiment is very satisfactory. In 
fact, the discrepancies for the excitation energies are all in the order of few tens of keV, the 
only exception being the J n = y + state, which comes about 170 keV below its experimental 
counterpart. Excluding the \ + state and the two levels at 1.35 and 2.65 MeV, for which we 
have not attempted any identification, the a value is only 64 keV. 

It should be noted that all the ten levels arising from the /z|/ 2 configuration lie at an 
excitation energy smaller than 1.5 MeV. In this energy region we also find the two seniority 
v = 1 states of the ^9/2/7/2 and ^9/2^13/2 configurations. The latter states as well as the 
ground state contain, however, significant configuration mixing. In fact, the percentage of 
configurations other than the dominant one is 22% in the ground and (| )i states, reducing 
to 14% in the (y + )i state. In all other levels up to 1.5 MeV the percentage of the dominant 
configuration ranges from 90 to 98%. Above 1.5 MeV the negative-parity states are mem- 
bers of the multiplet ^9/2/7/2, while all the positive-parity ones originate from the /ig/ 2 ^i3/2 
configuration, except the (y + )2 state, which arises from the ^9/2/7/2^13/2 configuration. All 
these states are essentially pure, the only exception being the (| )^ state, which contains 
38% of the ^9/2/7/2 configuration. 



C. Spectrum of 212 Rn 



Rather little experimental information [|T1| is presently available for 212 Rn. Up to about 
4 MeV only 22 excited states have been identified (nine of them with unknown spin and 
parity), while our calculations predict a much higher level density. In particular, in the 
low-energy region (up to 2.6 MeV) we find 27 states compared to 8 in the experimental 
spectrum. In this situation, any attempt to associate calculated states with experimental 
levels without assigned spin and parity may be misleading. Therefore, in Fig. 3 we exclude 
such states in the experimental spectrum and report only those yrast and yrare calculated 
states which are candidates for the observed levels. For completeness, all the calculated 
excitation energies up to about 2.6 MeV are listed in Table IV. It should be mentioned 
that above 4 MeV excitation energy several high-spin states have been observed. In Fig. 
3, however, we do not include these levels, since their description is likely to require that 
core-excited states be explicitly taken into account. 

From Fig. 3 we see that the calculated spectrum reproduces very well the experimental 
one, the discrepancies being smaller than 100 keV for the energies of 9 out of the 13 states 
considered. The rms deviation a is only 85 keV, in line with the values obtained for the two 
lighter isotones. 

As regards the structure of the states, we find that the wave functions of the seven 
higher-lying levels are substantially pure. These states are members of the three multiplets 
hl /2 (r = 10+, 12+), h 3 g/2 f 7/2 (J w = 12+14+), and h*i 13/2 (J* = 15" 16" 17"), and the 
percentage of the dominant configuration is at least 95%. This is not the case for the lower- 
lying states, whose wave functions contain significant configuration mixing. In the first four 
excited states the percentage of the dominant configuration, hg^ 2 , ranges from 71 to 76% 
while it becomes 55% in the ground state. As for the 8+ and 11~ states, the percentages of 
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the dominant configurations, ^9/2/7/2 and h^ 2 i 13 / 2 , are 80 and 82%, respectively. 



D. Electromagnetic properties 

The effective operators needed for the calculation of electromagnetic observables have 
been derived as described in Sec. II. 

In Table V we compare the experimental magnetic moments in 210 Po, 211 At and 212 Rn 
T8| , |36[| with the calculated values. We see that the agreement is remarkably good in all 



cases. Only two Ml reduced transition probabilities are known in 211 At |18|j37| . They are 
compared with our theoretical results in Table VII. We see that both the calculated and 
experimental values are extremely small. 

Let us now come to the electric observables. In Tables VI and VII we compare the 
calculated quadrupole moments and E2, E3 transition rates with the experimental ones 
r8] . pB| -p5| . The agreement is very good, the only discrepancy regarding the B(E2; 2f — > 0^) 



in 210 Po. It should be mentioned, however, that the experimental value was obtained by 
comparing the cross section cr(2 + , 210 Po) measured in a 210 Po study by inelastic scattering 
of deuterons with the corresponding one for 206 Pb f40 |. As far as the quadrupole moments 
are concerned, 4 out of the 5 calculated values are within the error bars and the observed 
signs, when measured, are correctly reproduced. It is worth noting that our results do not 
differ significantly from those obtained using an effective proton charge = 1.5e, which is 
consistent with the values adopted by other authors fI~6| , |T7| . 



As regards the _B(.E3)'s, they are all underestimated by our calculations. It is well known 
that enhanced E3 transitions in nuclei in the lead region can be taken as a signature of 



mixing of the 3~ core excitation into the involved levels ||41|| . We should note, however, that 
whereas our calculations fail to reproduce the B(E3) values, a good description of the states 
involved in such transitions is obtained for the excitation energies as well as for the other 
electromagnetic properties. Thus, these states are likely to contain very small components 
of octupole excitation which, however, are sufficient to largely enhance the E3 transition 
rates. In particular, the E3 transitions in 211 At and 212 Rn, and the 11^ — > 8^ transition in 
210 Po correspond to the single-particle transition z 13 / 2 — > fan, which is expected to be very 
fast owing to the coupling between the orbital and the 3~ collective state []4~2]| . The 
Hi ~ * &t transition in 210 Po is instead of the type — > ^9/2 and is slowed down by spin 
flip. 



IV. DISCUSSION AND CONCLUSIONS 

In this work, we have performed shell-model calculations for the N = 126 isotones 210 Po, 
211 At, and 212 Rn, employing an effective interaction derived from the Bonn A nucleon- 
nucleon potential by means of a G-matrix folded-diagram method. As for the single-proton 
energies, we have taken them from the experimental spectrum of 209 Bi. It should be stressed 
that, since we have also derived in a microscopic way the effective operators needed for 
the calculation of electromagnetic observables, no use has been made of any adjustable 
parameter. 
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These calculations, as well as the previous ones on neutron hole Pb isotopes fL4,21 
are the first in the 208 Pb region where a modern realistic interaction has been used. As 
already mentioned in the Introduction, the first attempt to employ in this region effective 
interactions derived from the free nucleon-nucleon potential dates back to the early 70s 
|T7|| . In that work, however, the Hamada- Johnston NN potential was used and only the 
3p-lh core-polarization diagram (the so-called bubble) was included in the calculation of the 
effective interaction. It should also be noted that to obtain good agreement with experiment 
for the 204_206 Pb isotopes, the bubble was multiplied by the empirical factor 0.75. 

As regards our calculations, we have obtained a very good description of both N = 126 
isotones and Pb isotopes in a truly microscopic way. It cannot be said, however, that our 



agreement with experiment is much better than that obtained in Ref . |17j • In this connection, 
we found it worthwhile to calculate the complete energy spectrum of 210 Po up to 5 MeV 
making use of the KH effective interaction. It turned out that the a value relative to the 
37 levels considered in Sec. Ill A is 116 keV, namely only about 30 keV larger than our 
value. On the other hand, it should be mentioned that a comparison between the results of 
the two calculations evidences more substantial differences. More precisely, the calculation 
with the KH interaction predicts some levels to lie 300-400 keV above those obtained with 
our effective interaction. We do not feel, however, that a detailed comparison between the 
two kinds of calculations is very meaningful. We consider as the main achievement of our 
studies of nuclei around 208 Pb to have shown that effective interactions derived from the 
Bonn A potential by means of a G-matrix folded-diagram approach lead to a quite accurate 
description of these nuclei. This outcome acquires more relevance when considered along 
with the results of our studies on nuclei with few valence particles or holes in the region 
of doubly magic 100 Sn and 132 Sn |T]-0j5|]. In fact, the remarkable overall agreement with 
experiment obtained in all cases considered leads to the conclusion that the new generation 
of realistic effective interactions is quite adequate for nuclear structure calculations. 

Actually, being focused on identical particle systems, our studies provide a stringent test 
of the isospin T = 1 matrix elements of the effective interaction. A test of the T = matrix 
elements is of course equally important. In this connection, it may be mentioned that in 



earlier works |12] it turned out that not enough attraction was provided by the calculated 
matrix elements of the T = effective interaction, which has a stronger dependence on the 
tensor force strength than the T = 1 interaction. We should point out, however, that in a 
recent study [f| of the doubly odd nucleus 132 Sb we have obtained results which are as good 
as those regarding like nucleon systems. Along the same lines we are currently studying 
other nuclei with both neutrons and protons outside closed shells. 

A main question relevant to microscopic nuclear structure calculations is the extent 
to which they depend on the NN potential used as input. We are currently trying to 
explore this problem. Preliminary calculations indicate that different NN potentials produce 
somewhat different nuclear structure results |]I|,|4"3"fl. In particular, it has turned out that the 
best agreement with experiment is produced by the Bonn A potential. 

In conclusion, at the present stage of our investigation of the role of realistic effective 
interactions in complex nuclei, it is our belief that a truly microscopic description of nuclear 
structure properties is now within reach. 
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FIGURES 

FIG. 1. Experimental and calculated spectrum of 210 Po. 
FIG. 2. Experimental and calculated spectrum of 211 At. 
FIG. 3. Experimental and calculated spectrum of 212 Rn. 
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TABLES 



TABLE I. Experimental and calculated ground-state binding energies (MeV) relative to 208 Pb 
for 210 Po, 211 At, and 212 Rn. 



Nucleus 


Expt. 


Binding energy 


Calc. 


210p o 


8.782 ± 0.004 




8.789 


211 At 


11.765 ± 0.005 




11.816 


212 Rn 


16.065 ± 0.006 




16.146 
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TABLE II. Comparison of the experimentally observed spectroscopic strengths from the 
209 Bi(a, i) 210 Po and 209 Bi( 3 He,<i) 210 Po reactions with the calculated values. See text for com- 
ments. 



l j 


J 


-^exp(.J-VleVJ 


^calc(,J-VleVJ 


(2J+1) a(n ^ 

(2,j l +i) y a ^ 1 ) 


(2J+1) cv3u„ a\ 

(2j7+i) b ( He ' d ) 


uaic. 


7 

"9/2 


2^ 


1 1 O 1 

1.181 


i i on 

1.130 


1.00 ± 0.06 


i i /' in nn 

1.16 ± 0.22 


0.98 


A + 

4^ 


1.427 


1.395 


i i n n^r 

1.79 ± 0.07 


i r o i n no 

1.58 ± 0.28 


1.78 






1.473 


1.493 


2.64 ± 0.09 


n o i n or 

2.63 ± 0.25 


2.58 




8^ 


1.557 


1 CT V CT 

1.555 


o A n i n ~i i 

3.40 ± 0.11 


3.42 ± 0.26 


O ■ ) (■ • 

3.36 


/V/2 


8 + 


2.188 


2.179 


1.71 ± 0.03 


1.64 ± 0.05 


1.68 


2+ 


2.290 


2.292 


0.42 ± 0.01 


0.42 ± 0.02 


0.47 




6+ 


2.326 


2.367 


1.31 ± 0.03 


1.26 ± 0.04 


1.28 




4+ 


2.382 


2.394 


(0.90) 


(0.90) 


0.88 




1 + 


2.393 


2.220 


(0.31) 


(0.35) 


0.30 




5 + 


2.403 


2.422 


/ 1 1 n\ 

(1.10) 


(1.10) 


1.10 




3 + 


2.414 


2.380 


(0.72) 


(0.75) 


0.70 




7^ 


2.438 


2.437 


i n i n no 

1.51 ± 0.03 


i rn i n r\r* 

1.50 ± 0.06 


i f n 

1.50 


H3/2 


3~ 


2.846 


2.862 






0.69 




2.849 


2.700 


3.25 ± 0.06 




2.30 




5 


2.910 




n o 1 i n m 

0.31 ± 0.01 








9~ 


2.999 


3.016 


(1.88) 




1.89 




7~ 


3.016 


3.065 


(1.53) 




1.50 




2~ 


3.024 


2.682 


(0.54) 




0.50 




5~ 


3.026 


3.024 


(0.78) 




1.10 




4r 


3.075 


3.039 


0.78 ± 0.02 




0.90 




6~ 


3.125 


3.097 


(1.34) 




1.30 




8" 


3.168 


3.121 


(1.66) 




1.70 




10 _ 


3.183 


3.154 


2.11 ± 0.04 




2.10 


£ 

J 5/2 


2 T 


3.792 


3.828 


n or i n no 

0.35 ± 0.02 




0.34 


A -I- 

4+ 


4.027 


4.152 


0.60 




0.67 




6 + 


4.139 


4.256 


0.82 




0.86 




n-4- 

3+ 


4.320 


4.309 


0.86 ± 0.04 




0.69 




5 + 


4.382 


4.391 


1.19 ± 0.05 




1.10 




6 + 


4.469 


4.503 


0.55 ± 0.03 




0.43 




7+,4+ 


4.553 


4.384,4.552 


1.83 ±0.07 




1.50,0.20 


P3/2 + Pl/2 


4+ 


4.027 


4.152 




0.03 


0.07 




6+ 


4.139 


4.256 




0.20 


0.29 




6+ 


4.469 


4.503 




0.56 ±0.04 


0.66 




4+ 


4.553 


4.522 




0.35 ±0.07 


0.22 




3+ 


4.591 


4.605 




(0.75) 


0.69 




5+ 


4.624 


4.673 




(1.35) 


1.10 




6+ 


4.644 


4.691 




(0.55) 


0.31 
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TABLE III. Calculated low-energy levels of 211 At. 





E(MeV) 


J 71 " 


/TV /T T T\ 

E(MeV) 


J w 


T" 1 /tv /r t r\ 

E(MeV) 


9 — 
2 


0.0 


5 — 
2 


2.040 


19 — 
T 


2.350 


7 — 
2 


0.679 


7 — 
2 


2.042 


9 + 
2 


2.393 


7~ 
2 


0.783 


11 — 

T 


2.042 


11 + 

T 


2.412 


5 — 
2 


0.955 


17 — 
~2~ 


2.045 


19 + 
TT 


2.427 


13 — 
~2~ 


1.053 


9 — 
2 


2.080 


29 + 
T 


2.466 


11 — 
~2~ 


1.098 


11 — 
T 


2.110 


21 + 

T 


2.472 


3 ~ 
2 


1.103 


l — 

2 


2.110 


25 + 

T 


2.515 


9 — 
2 

13 + 
~2~ 


1.186 


15 + 

2 


2.124 


23 + 
T 


2.528 


1.236 


3 — 
2 


2.131 


11 + 
T 


2.530 


15 — 
T 


1.337 


13 + 
T 


2.136 


7 + 
2 


2.539 


17 — 
~2~ 


1.339 


21 — 
~2~ 


2.178 


3 + 
2 


2.541 


21 — 
~2~ 


1.467 


7- 
2 


2.181 


5 + 
2 


2.549 


9 — 
2 


1.631 


5 — 
2 


2.185 


11 — 

T 


2.560 


7 — 
2 


1.681 


9 — 
2 


2.189 


9 + 
2 


2.574 


11 — 
~2~ 


1.721 


13 — 
T 


2.204 


13 + 
2 


2.589 


3 — 
2 


1.824 


9 — 

2 


2.229 


27 + 
2 


2.609 


5- 
2 


1.856 


11- 
2 


2.238 


5- 
2 


2.617 


13- 
2 


1.929 


15- 
2 


2.245 


13- 
2 


2.625 


15- 
2 


1.940 


17 + 

2 


2.279 


17 + 
2 


2.629 


9- 
2 


1.967 


13- 
2 


2.282 


9- 
2 


2.686 


23~ 
2 


1.969 


17- 

2 


2.289 


9 + 
2 


2.688 


19- 
2 


1.987 


15- 
2 


2.291 
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TABLE IV. Calculated low-energy 


l i c 212r> 

levels oi Rn. 




r 


E(MeV) 


J* 


E(MeV) 


0+ 


0.0 


3+ 


2.276 


2+ 


1.221 


4+ 


2.277 


4+ 


1.506 


5+ 


2.290 


6+ 


1.619 


5+ 


2.357 


8+ 


1.677 


7+ 


2.387 


4+ 


2.057 


4+ 


2.450 


8+ 


2.122 


3+ 


2.473 


0+ 


2.170 


0+ 


2.499 


6+ 


2.177 


6+ 


2.561 


2+ 


2.198 


2+ 


2.581 


1+ 


2.208 


11" 


2.597 


2+ 


2.211 


8+ 


2.632 


6+ 


2.226 


3" 


2.651 


7+ 


2.274 


10+ 


2.655 
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TABLE V. Calculated and 


experimental dipole moments (m nmj. 




Nucleus 


r 


A* 

Calc. 


Expt. 


210p o 


6f 


+5.29 


±5.48 + 0.05 






+7.06 


+7.35 + 0.05 




nr 


+ 13.12 


+12.20 + 0.09 


211 At 


(F)i 


+6.6 


±6.8 + 0.6 






+9.32 


+9.56 ± 0.09 




(f + )l 


+ 16.23 


+15.31 + 0.13 


212 Rn 


4f 


+3.56 


±4.04 ± 0.24 




6f 


+5.308 


±5.454 ± 0.048 




8f 


+7.064 


+7.152 + 0.016 




14+ 


+15.07 


±14.98 ± 0.42 




i7r 


+ 18.45 


±17.85 ±0.17 



TABLE VI. Calculated and experimental quadrupole moments (in eb). 



Nucleus 


r 


Calc. 


Q 


Expt. 


210p o 




-0.588 




-0.552 ± 0.020 




ur 


-0.92 




-0.86 + 0.11 


211 At 


(f")i 


-0.54 




±0.53 ± 0.05 




(f + )i 


-1.07 




±1.01 ±0.19 


212 Rn 


8 + 


-0.29 




(-)0.17±0.02 
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TABLE VII. Calculated and experimental reduced transition probabilities (in W.u.). 



Nucleus 


Transition 


TTT , TTT 

J i -» J f 


Reduced transition probabilities (in W.u.) 






Calc. 


Expt. 


210p o 


E2 


2 + -0 + 


3.55 


0.56 ±0.12 




E2 




4.46 


4.53 ±0.15 




E2 


6 + -4+ 


3.07 


3.00 ±0.12 




E2 


8 + -6+ 


1.25 


1.10 ±0.05 




E3 


Hi ^8+ 


6.1 


19.7 ± 1.1 




E3 




0.55 


3.71 ±0.10 




E2 


/ 3 — \ / 5 — \ 
(§ )l^(ll) 


10.1 


12.5 ± 1.4 




E2 


/ 3 ~ \ c 7 — \ 


1.67 


1.77 ±0.17 




E2 


(1 )l 


0.15 


0.39 ± 0.04 




E2 


(¥ )i-(t )i 


2.3 


1.3 ±0.3 




E2 


(t )i-(¥ )i 


2.60 


2.51 ±0.05 




E2 


(f + )l-(f + )l 


1.6 


1.8 ±0.5 




Ml 


(r)i-d _ )i 


8 x 10~ 7 


1.4 x 10~ 4 ±0.4 x 10~ 4 




Ml 


(F)i-(F)i 


8 x 10~ 7 


0.7 x 10" 4 ±0.2 x 10~ 4 




E3 


(f + )i-(r)i 


6.3 


20.1 ± 1.8 


212 Rn 


E2 




1.42 


1.04 ± 0.04 




E2 


6+^4+ 


0.73 


0.40 ± 0.05 




E2 




0.252 


0.115 ±0.006 




E2 


12+ -» 10+ 


3.6 


4.4 ±0.2 




E2 


14+ -» 12+ 


0.008 


0.032 ± 0.008 




E2 


14+ -» 12+ 


3.4 


3.6 ±0.5 




E2 


i7r - isr 


2.9 


3.0 ± 1.6 




E3 


177 - 14+ 


6 


16 ±6 



20 



211 



At 



(27/2+) 



27/2+ 



> 

CP 



,(27/2+) 
■1/2,3/2 
■(29/2) + 
■(25/2+) 
.(1/2) + 



2 - 



1 - 



.-(5/2-) 
-(3/2)- 
'-■1/2 .3/2 
-(5/2)- 
'■■■(23/2>- 
-(3/2)- 

,(21/2)- 
',(9/2,11/2,13/2) 

:.-(i7/2)- 

■ ■■(15/2 ,17/2 ) 
,(11/2,13/2)- 
-(3/2)- 
■■(13/2)- 
-(5/2)- 
-(7/2)- 

(7/2)" 



,3/2+ 

,3/2+ 

,1/2+ 

,27/2+ 

■-3/2+, 

,25/2+ 

,29/2+ 

:5/2- 

■;3/2 

,172- 

'.■7/2 _ 

',5/2" 

,19/2- 

,23/2 _ 

-9/2- 

-15/2- 

■13/2" 

■5/2- 

■3/2-_ 

■■11/2 

-7/2- 

■9/2- 

■ 21/2" 
-17/2- 
■15/2- 
-13/2+ 

■9/21 
-3/2 

■ ll/2- 
■13/2- 
■5/2- 

-7/2" 
-7/2" 



- 



9/2" 



9/2" 



Expt. 



Calc. 



